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Abstract 

The results from a concentrated computational ef- 
fort are presented with the primary objective being di- 
rected at evaluating the vortical-flow-prediction capa- 
bility of an unstructured-grid Euler solver. Both vis- 
cous and inviscid solutions, obtained from an estab- 
lished structured-grid method, along with an experi- 
mental wind-tunnel data are used as bench-mark mear 
sures to assess the validity of the unstructured-grid Eu- 
ler results. Viscous effects on vortical flows are first 
identified by comparing the viscous and inviscid solu- 
tions obtained from the structured-grid method. Com- 
putational data analysis are then presented which re- 
veal excellent correlations between the inviscid struc- 
tured and unstructured-grid results in terms of off- 
surface flow structures, surface pressure distribution 
and the predicted longitudinal aerodynamic character- 
istics. The sensitivity of the unstructured-grid inviscid 
solutions to grid refinement is also discussed along with 
an analysis of the convergence and performance char- 
acteristics for each method. 


Nomenclature 

Cd drag coefficient, DTa.g/q 00 S re f 

C Do measured drag coefficient at Cl = 0 
Cl lift coefficient, Lift/7 co 5 r «/ 

C m pitching moment coefficient, 

pitching moment/5oo5 re /c re /, 
moment reference center at 1.86c re y 
C p pressure coefficient, (p-Poo)/?oo 

c re f reference chord, 10.92 inches 

FS fuselage station 

Moo freestream Mach number 

MTVI modular transonic vortex interaction 
p static pressure 

p 0 free-stream total pressure 

Po t i local total pressure 

Po,}/Po normalized total pressure 

Poo freestream static pressure 

qoo freestream dynamic pressure 

R/t Reynolds number per foot 

Sref area of reference wing planform, 104.1 in 2 

v* wall-friction velocity, \Jr w j p 

X longitudinal fuselage station, inches 
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y + inner-law variable, yv* fv 

a angle of attack, degrees 

p viscosity 

u kinematic viscosity, p/p 

p density 

r w wall shear stress 

Introduction 

The unstructured-grid approach is an emerging 
new technology with a tremendous potential to ease the 
application of Computational Fluid Dynamics (CFD) 
methodologies for solving the flow about complex con- 
figurations in a timely manner. Unlike the conven- 
tional structured-grid approach which uses hexahedral 
elements to discretize a computational domain about a 
given aerospace vehicle, this new technology utilizes un- 
structured tetrahedral meshes. The inherent flexibility 
with which a tetrahedral mesh can conform to a ge- 
ometrically complex configuration is generally consid- 
ered to be the primary advantage of the unstructured- 
grid approach over the structured-grid method. How- 
ever, as the initial development of the technology is 
maturing, other advantages have been identified such as 
solution-adaptive grids. Progress has been made in the 
past few years in both the tetrahedral-mesh generation 
techniques (Refs. 1-4) as well as solution algorithms 
(Refs. 5-9). However, unstructured-grid methodology 
still has some fundamental shortcomings which hinder 
its full potential applicability to complex configuration 
across a wide range of flow conditions. 

One of the major impediments to a wide spread 
application of this new technology is presently due to 
its restricted viscous-flow capability. Some progress has 
been made in computing the viscous flow about 2-D air- 
foils and with limited success for simple 3-D configura- 
tions. These viscous flow calculations are generally per- 
formed on so called hybrid grids which incorporate cer- 
tain structured-mesh arrangements near the surface to 
facilitate the boundary-layer flow-resolution. An exam- 
ple of such grid is the prismatic/tetrahedral mesh (Refs. 
10-11) arrangement where prism cells are generated, in 
a quasi-structured fashion, from a triangulated surface 
grid definition. These prisms are generally confined to 
the viscous layer region beyond which a conventional 
unstructured grid generation method such as advanc- 
ing front (Refs. 1, 2) or Delaunay triangulation (Refs. 
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12, 13) can be utilized to fill in the inviscid field without 
constraints. Another example of a hybrid grid is to sim- 
ply use a structured grid to resolve the boundary-layer 
flow and couple that to an inviscid unstructured-grid 
(Refs. 14, 15). In general, it appears that the obvi- 
ous challenge for extending the present unstructured 
grid technology to simulate viscous flow is dependent 
on the development of an innovative approach for gen- 
erating highly stretched meshes. It would also be ad- 
vantageous for any new unstructured-grid generation 
approach to include directional features to accommo- 
date some of the existing algebraic turbulence models 
that are presently in use in structured-grid methodolo- 
gies. 

Recently there have been a number of success- 
ful inviscid-flow analyses performed about fairly com- 
plex configurations using the existing state-of-the-art 
unstructured-grid technology (Refs. 16, 17). These ap- 
plications have clearly demonstrated the inherent flexi- 
bility associated with tetrahedral meshes to conform, as 
well as being generated in a timely manner, about com- 
plex configurations. These advantages coupled with the 
inviscid-flow approximation (requires less intense com- 
putations, in terms CPU-time, compared to viscous 
flow simulations) provide an efficient computational 
method that can be used for the flow analysis of certain 
aerodynamic problems for which the viscous effects are 
known to be small, such as vortex-flow separation from 
a sharp-edged configuration. The latter computational 
method can also be used in the early configurational 
design to study the relative aerodynamic performance 
and trends due to a change in the geometrical compo- 
nents. Although, the inviscid computational-method 
provides a viable engineering tool for the preliminary 
flow analysis of a given aerodynamic problem, partic- 
ular attention should be paid for the interpretation of 
the final results and their dependencies on the fluid vis- 
cosity. As a result, the viscous effects have to be clearly 
quantified in order to eliminate the element of uncer- 
tainty which is always associated with an inviscid-flow 
analysis, in particular when correlated with experimen- 
tal data. 

The objective of the present investigation is to 
address the validity of the inviscid results obtained, 
from an algorithm that utilizes tetrahedral-mesh to dis- 
cretize the computational domain, at conditions that 
are conducive to flow separation leading to vortical- 
flow structures. The present study takes advantage of 
the only possible approach which is to use the com- 
putational results from an established structured-grid 
methodology as a bench mark to assess the validity 
of the unstructured-grid solutions. The present code- 
calibration strategy is broken into three steps as fol- 
lows; (1)- Assess the accuracy of the structured-grid 


results by comparing the turbulent, thin-layer Navier- 
Stokes solutions with experimental data. (2)- Isolate 
the effect of viscosity by complementing the above vis- 
cous computations with an inviscid analysis within the 
structured-grid framework. (3)- Assess the validity 
of the inviscid unstructured-grid results by compar- 
ing them to those obtained with inviscid structured- 
grid. Both the structured- and the unstructured-grid 
flow-solvers used in the present investigation, known 
as CFL3D (Refs. 18-21) and USM3D (Ref. 9) respec- 
tively, are developed at NASA Langley Research Cen- 
ter. 

Configuration, of Interest 
The configuration used in the present code cali- 
bration study is the isolated fuselage component of a 
wind-tunnel model, referred to as Modular Transonic 
Vortex Interaction (MTVI) model, shown in figure 1(a). 
The complete MTVI model employs a 60° sharp-edged 
cropped- delta wing with segmented leading-edge flap, 
and a chine shaped fuselage. In addition, the model is 
built to be fitted with two different vertical tail arrange- 
ments; a single centerline tail or an alternative twin 
outboard tails as shown in the figure. The model has 
been tested at various flow conditions, to investigate 
the interactions between the chine-forebody vortices 
and different vertical tail arrangements, in the 7- by 
10-foot High-Speed-Tunnel at NASA Langley Research 
Center. In addition to the wind-tunnel testing of the 
full MTVI configuration, several runs were made with 
the isolated fuselage primarily to support this compu- 
tational effort. A photograph of the isolated MTVI 
fuselage wind-tunnel model, mounted on the sting, is 
shown in figure 1(b). Two typical flow condition, listed 
in table 1, were selected from the latter test for the 
present code- calibration study. 


a 


R/t x 10 6 

19.8° 

0.4 

2.4 

29.8° 

0.4 

2.3 


Table 1. Selected flow-conditions. 


The rational for simplifying the MTVI configura- 
tion for the present code calibration study is primar- 
ily dictated by the geometrical complexity associated 
with the complete MTVI configuration and the foreseen 
difficulties with the structured-grid application part. 
These difficulties include the development of a multi- 
block grid-strategy that is suitable to the existing flow 
solver as well as the actual flowfield grid generation. 
These difficulties actually prompted the present code- 
calibration study to assess the vortical-flow-prediction 
capability of the unstructured-grid methodology so that 
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it could subsequently be applied to the full MTVI con- 
figuration with confidence. The isolated MTVI fuse- 
lage geometry incorporates sharp leading-edges on both 
the chine forebody as well as the exposed wing-fuselage 
juncture. Such surface discontinuities, introduced by 
the sharp leading-edges, are the principle cause of flow 
separation which subsequently leads to the formation of 
the vortical-flow structure. This type of flow separation 
is generally known to be insensitive to viscous effects 
which make it particularly appropriate for the present 
inviscid code- calibration study of vortical flows. 

Experiment 

The wind-tunnel experiment is conducted in the 
NASA Langley 7- by 10-foot High Speed Tunnel (Refs. 
22, 23). This is a closed circuit, continuous-flow at- 
mospheric tunnel with a solid wall test section 6.6 feet 
high, 9.6 feet wide, and 10 feet long. The tunnel oper- 
ational Mach number range from 0.0 to 0.9 with maxi- 
mum Reynolds number of about 4 million per foot. 

The MTVI wind-tunnel model is instrumented to 
measure surface static pressures at three longitudinal 
stations on the forebody and three longitudinal sta- 
tions on the aft fuselage that are extended outboard 
onto the wing surface. The pressures are measured on 
both sides (i.e., starboard and port) of the model de- 
signed to address the flow asymmetry, if any. These 
pressure stations are highlighted in Fig. 1(a). The 
sting-mounted wind-tunnel model was equipped with 
an internally mounted strain gauge balance to measure 
the six-component forces and moments. The wind- 
tunnel testing of the isolated MTVI fuselage (see fig- 
1(b)) uses similar mounting apparatus for measuring 
the forces and moments. However, surface static pres- 
sures are measured only at the three longitudinal sta- 
tions on the forebody. These fuselage stations (FS) are 
located at X=6.1, 10.45, and 14.5 inches measured from 
the nose apex. For reference, the total fuselage body 
length is 32.48 inches. 

The measured wind-tunnel data are corrected for 
effects associated with angle-of-attack, wall interference 
and the model base pressure. The model support sys- 
tem incorporated a pitch sensing device to measure the 
angle-of-attack which is subsequently corrected to ac- 
count for the balance and sting deflection under load. 
The wall interference effects are accounted for through 
the notions of blockage (ref. 24) and jet boundary (ref. 
25) correction principles. The model-base pressures are 
measured for correcting the data to exclude the result- 
ing force acting on the base of the model. The latter 
correction is employed to re-adjust the base pressure to 
that of the freestream static pressure. 

Computational Methods 

This section provides a brief background on the 


computational methods used in the present investiga- 
tion along with some of their respective capabilities. 

CFL3D: The unsteady Reynolds- averaged Navier 
Stokes equations are solved in a body-fitted coordi- 
nate system. The algorithm is based on cell centered, 
finite volume approach incorporating upwind-biased, 
flux- difference splitting for the inviscid terms. Sev- 
eral turbulence modeling capabilities are implemented, 
however, the present computational effort only uses the 
Baldwin k Lomax model (Ref. 26) with Degani k 
Schiff (Ref. 27) modification to account for massive 
crossflow separation. The code capabilities and some 
of the new improvements made over the past few years 
are evaluated and discussed in references 18-21. For 
applications to complex configuration, the method has 
been recently extended to multiblock structured grids 
with the concept of generalized surface patching (Ref. 
21) which has been successfully applied to the F/A-18 
aircraft configuration (Ref. 28). 

USM3D: The time dependent, three dimensional Eu- 
ler equations are solved in a computational domain 
discretized by tetrahedral-mesh elements. Similar to 
CFL3D, USM3D is based on cell centered, finite vol- 
ume approach which also uses an upwind flux- difference 
splitting method. An implicit Gauss-Siedel scheme 
has recently been implemented to advance the solution 
(Ref. 29). Convergence to steady state is accelerated 
with a local time stepping. Further details are discussed 
in reference 9 and the results from some recent appli- 
cations can be found in references 16 and 17. 

Computational Grid Attributes 

Surface Definition: The isolated MTVI fuselage con- 
figuration is 32.48 inches long, 4.80 inches wide and 3.00 
inches thick. An oblique view of one-half the configura- 
tion surface definition used in both viscous and inviscid 
structured-grid analysis is shown in figure 2. The com- 
plete geometry is defined analytically with the excep- 
tion of a small area, referred to as the “blended region’* 
shown with dotted lines in the figure, which will be dis- 
cussed later. The forebody portion is 15.08 inches long 
and defined in two parts with a longitudinal dividing 
cross-section at X=10.45 inches, where X is measured 
from the nose apex. The planform and the profile for 
the front part of the forebody is defined analytically by 
a circular-arc while the corresponding cross-sectional 
geometry is defined by a cubic relationship. The geom- 
etry for the back part of the forebody is defined by lon- 
gitudinally repeating the base cross-section of the front 
part up to where the wing leading-edge intersects the 
fuselage (i.e., X=15.08). The forebody cross-sectional 
definitions are spaced longitudinally so that there are 
three stations that match those of wind-tunnel model 
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pressure stations. These fuselage stations are denoted 
in fig. 2. 

The fuselage geometry aft of the forebody is gener- 
ated by modifying the chine leading-edge definitions to 
incorporate the streamwise thickness associated with 
the wing-fuselage juncture. The forward part of the 
wing-fuselage juncture is defined by a circular-arc up to 
where the wing maximum thickness occurs (X=21.08), 
beyond which it remains constant. The fuselage geom- 
etry aft of the forebody was then completed by adding 
the surface definition for the blended region. That in- 
formation was extracted from the Computer Aided De- 
sign (CAD) surface patch definition used to build the 
wind-tunnel model. This blended region describes the 
area where the leading-edge thickness, (Le., added to 
the chine sharp- leading-edges to accommodate for the 
wing thickness) is smoothly blended inboard with the 
analytically defined cross-sections. 

The computational surface grid definition for the 
isolated MTVI fuselage is defined with two networks; 
the first includes the forebody geometry whereas the 
second consists of the aft fuselage. Each network has 
35 longitudinal stations. There are 61 and 69 circum- 
ferential grid points per station in the first and second 
network, respectively. A total of 4,550 grid points are 
used to define the computational surface geometry. The 
grids are clustered, with a cosine-like spacing, near all 
the sharp leading-edges, plane of symmetry and the 
network boundaries. Furthermore, to ensure smooth 
outflow conditions at the base, a third network was 
added which is nominally a body-length downstream 
extension of the configuration base cross-section defined 
with 18 longitudinal stations. 

Structured Flowfield-Grid: The body slenderness 
associated with the MTVI fuselage geometry is pri- 
marily used to guide the fiowfield grid-strategy and 
topology selection. This led to the selection of a 
longitudinally-patched multiblock grid-strategy, with 
a C-0 topology for the forebody block, an H-0 for 
the aft fuselage and the downstream grid extension. 
The fiowfield grid is generated with an in-house de- 
veloped code that is based on transfinite interpolation 
techniques and have been successfully applied to vari- 
ous configurations in the past (Refs. 30, 31). A tur- 
bulent fiowfield grid is first generated with sufficient 
normal clustering near the surface to adequately re- 
solve the laminar sublayer associated with the turbu- 
lent boundary layer flow at the subject freestream con- 
ditions (a = 19.8®, Moo = 0.4, R/t — 2.4 x 10 6 ). This 
grid produced an average normal cell size next to the 
wall of approximately 0.0005 inches (i.e., 1.5 x 10 -5 /, 
where l is the total configuration body length) which 
corresponds to y + ft 3. The farfield boundary is located 


approximately one body length away in the radial di- 
rection. 

The viscous computations, with turbulent-flow as- 
sumptions, are first performed in an effort to determine 
the approximate thickness of the boundary-layer. A de- 
tailed flow analysis of the converged turbulent solutions 
revealed that the boundary-layer flow is resolved with 
about 18 cells near the plane of symmetry at a mid- 
body station. Subsequently, the inviscid fiowfield grid 
is generated by redistributing the points radially within 
the existing viscous grid blocks. The inviscid nearfield 
grid is generated by ensuring that an average normal 
cell size, next to the wall, approximately matched the 
edge of the turbulent boundary-layer at a mid-body sta- 
tion. At this station, the inviscid normal cell size next 
to the wall is determined to be about 0.1 inches (i.e., 
3.1 x 10 -3 /) which can roughly occupy about 24 viscous 
cells. Figure 3, shows an isometric farfield (a), nearfield 
(b) views of both the viscous (dotted lines) and invis- 
cid (solid lines) structured-grids that are superimposed 
over one another. As evident in Fig. 3(b), the inviscid 
fiowfield grids are generated in such a manner that the 
nearfield grid density (region where the primary vor- 
tex structure is expected to occur) remain roughly the 
same as those of the viscous grid. The final viscous 
and inviscid fiowfield grids contained 65 and 49 points 
in the radial direction, respectively. Similarly, there are 
a total of 376,480 and 283,808 points in the complete 
viscous and inviscid fiowfield grid, respectively. 

Unstructured Flowfield-Grid: The strategy used to 
generate the unstructured grid about the MTVI fuse- 
lage configuration is primarily guided by the structured 
grid analysis. These strategies include the use of the 
existing structured surface grid definition as the data 
base for generating the surface triangles. Similarly, the 
fiowfield grid boundaries such as the plane of symme- 
try, farfield and the outflow are chosen to be identi- 
cal to those used in the structured grid analysis. The 
surface and the fiowfield grid boundaries are initially 
triangulated. The desired grid distribution as well as 
the spacing on both the surface and the fiowfield grid 
boundaries are controlled by a set of parameters deter- 
mined from coupling a uniform Cartesian background 
grid with a series of point/line source elements (Ref. 
32). The surface triangular elements are generated 
so that the cell density and spacing, such as around 
the leading-edges, are comparable to those generated 
with the structured grids. A similar approach is also 
adopted to generate the desired grid density and spac- 
ing in the fiowfield boundaries (i.e., plane of symmetry, 
farfield, and the outflow plane) comparable to the struc- 
tured grids. Subsequently, the triangulated surface and 
the fiowfield boundaries are used to generate the corre- 
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sponding volume tetrahedral-mesh with the advancing 
front method known as VGRID3D (Ref. 3). Figure 5 
presents the final unstructured grids from a similar van- 
tage point as the structured-grid arrangement shown 
earlier in figure 3. 

Despite much attention being paid to minimizing 
the differences between the unstructured and its coun- 
terpart structured grid, in terms of density and spac- 
ing, the final grids still contain some dissimilar char- 
acteristics. One of the distinct dissimilarities can be 
attributed to the cell aspect-ratio. Unlike structured 
grids, the unstructured grid generation techniques, such 
as the present advancing front, are primarily designed 
to generate cells having a unit aspect ratio. In the 
present investigation, this limitation coupled with the 
constraints in the desired grid density and spacing, was 
foreseen to lead into an unacceptable number of cells 
in the final volume grid, which would have about six 
times more cells than the structured-grid. In an ef- 
fort to alleviate this discrepancy, the present unstruc- 
tured grids are generated with an aspect ratio of about 
three. The effect of this simple grid stretching in the ax- 
ial direction was significant enough to reduce the total 
number of cells by a factor of about three. The final 
unstructured-grid domain consisted of about 102,000 
nodes and 567,000 cells which is roughly twice as many 
cells as used in the inviscid structured-grid. The pri- 
mary reason for the excess number of tetrahedra is re- 
lated to the differences in cell aspect-ratio. These dif- 
ferences can clearly be seen by comparing the surface 
grid stretching in the axial direction in Figs. 4 and 
2. The present investigation did not seek to generate a 
tetrahedral mesh with higher aspect-ratio due to 1)- the 
anticipated difficulties with grid generation and 2)- the 
sensitivity of solutions from USM3D to highly stretched 
tetrahedral meshes has not been fully addressed. 

Results and Discussion 

The results from the present code calibration study 
will be discussed in three categories with intention of 
highlighting the vortical flow prediction capabilities of 
the present unstructured-grid methodology. These cat- 
egories include the off-surface flow simulation, surface 
pressure and the force and moment predictions. In 
each category the inviscid computational results ob- 
tained from the unstructured-grid method will be corre- 
lated with complementary solutions acquired with the 
structured-grid method for both the inviscid as well 
as the thin-layer Navier-Stokes with turbulent flow as- 
sumptions. 

Off-surface flow correlations: A composite image 
summarizing the computational results in terms of nor- 
malized total-pressure contours ( Po,t/po)> obtained at 
a ~ 19.8°, Moo — 0.4, Rft = 2.4 x 10° from both 


the structured and unstructured-grid methodologies, 
are presented in figure 6 for various cross-flow planes. 
The figure has been constructed such that the verti- 
cal plane of symmetry, relative to the isometric view 
of the surface grid representation as shown in the fig- 
ure, divides the computational results into two parts; 
the structured-grid solution part (i.e., starboard) and 
the unstructured-grid solution part (i.e., port side). 
Coupled with the surface grid representations are the 
corresponding inviscid solutions (i.e., gray-scaled nor- 
malized total-pressures) shown in six cross-flow planes 
along the length of the isolated MTVI fuselage. For 
a more quantitative comparison, the figure also shows 
the predicted normalized total-pressure contours as well 
as the corresponding cross-sectional geometries, from 
a head-on vantage point, for the first three stations 
on the forebody. To assess the viscous effects on the 
resulting vortex-flow structures, the turbulent-flow so- 
lutions computed by the structured-grid method are 
also included for the same three cross-sectional geome- 
tries. It should be noted that these three forebody sta- 
tions are chosen to be identical to those experimen- 
tally measured pressure stations and also the normal- 
ized total-pressure contours are plotted over the same 
range (( Po,t/Po)min - 0.86, (p 0 ,i/p o )max ~ 1-0) and 
levels ((, Po,llPo)inc ~ 0*0175). 

The normalized total-pressure contours obtained 
from the structured-grid viscous solutions indicate a 
flow structure that can generally be characterized by 
three features. These flow features include the primary 
and secondary vortex flow structures as well as a thin 
boundary-layer region near the surface. The flow sepa- 
ration from the sharp leading edges of the chine config- 
uration is the principal cause for the formation of the 
primary vortex flow. The core (here defined as the min- 
imum normalized total-pressure) of the primary vortex 
appears to move up and slightly outboard with increas- 
ing axial distance. Also, note that the normalized total- 
pressure levels in the primary vortex core remain essen- 
tially the same with the exception of a slight increase 
at FS 3. As expected, the normalized total-pressure 
contours obtained from the structured-grid inviscid so- 
lutions reveal only the primary vortex flow structure 
(see figure 6). Relative to the viscous solutions, the 
predicted inviscid primary vortex appears to have a 
similar axial core-trajectory, but with a contour-level 
structure that is more compact and clustered. 

The corresponding normalized total-pressure con- 
tours obtained from the unstructured-grid inviscid solu- 
tions are shown in the right hand-side column of figure 
6 for the same three forebody cross-sectional geome- 
tries. As evident, the unstructured-grid inviscid predic- 
tions of the primary vortex-flow structure are incredi- 
bly close to those computed with the structured-grid 
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inviscid method in terms of axial core-trajectory and 
the total-pressure contour levels. This excellent cor- 
relation between the structured- and unstructured-grid 
inviscid solutions for predicting the primary vortex-flow 
also persists for the aft-fuselage stations. For the latter 
stations, the corresponding solutions (i.e,, gray-scaled 
normalized total-pressures), in cross-flow planes, are 
augmented with the isometric view of the surface-grid 
representations in figure 6. 

Surface pressure correlations: The computed span- 
wise surface pressure coefficients obtained from struc- 
tured grid methodologies, based on both turbulent and 
inviscid flow simulations, along with the unstructured- 
grid inviscid results are presented in figure 7 for the 
same three stations on the forebody as those used in the 
previous discussions. In addition, the figure presents 
the corresponding experimental data measured on both 
the starboard and port side of the model. Subsequent 
analysis of the data presented in figure 7 addresses: 1)- 
the accuracy of the structured- grid viscous solutions in 
predicting the experimental data; 2)- the viscous effects 
from the structured-grid solutions; and 3)- the validity 
of the unstructured-grid inviscid solutions. 

The surface pressure coefficients obtained from the 
structured-grid viscous solutions, with turbulent-flow 
assumptions, correlate reasonably well with the experi- 
mental data for both the upper and lower surface. The 
agreement appears to be particularly good at FS 2 and 
FS 3 in the suction-peak region associated with the pri- 
mary vortex flow. There is also a good surface pressure 
correlation, at FS 3, in the region where the secondary 
vortex-flow occurs (located just inboard of the lead- 
ing edge, see figure 6). However, this agreement is not 
quite as good at FS 2. The largest surface pressure 
disagreements between the computational results and 
data appear to occur at FS 1, where the primary vor- 
tex suction peak is underpredicted. The source of this 
disagreement has not yet been identified, but it could 
be associated with transitional flow in the experiment 
which is not computationally modeled. 

The viscous effects on the surface pressure co- 
efficients are clearly evident, figure 7, by correlating 
the two computational solutions (i.e., solid and dashed 
lines) obtained from the structured-grid method. The 
correlations indicate that the inviscid solutions develop 
a primary vortex suction peak that is higher (more neg- 
ative) and located slightly outboard compared to the 
viscous results. Unlike the viscous solutions, the com- 
puted inviscid surface pressure distribution do not in- 
dicate the presence of an inflection point that leads to 
a region of mild expansion which is directly attributed 
to the secondary vortex-flow effects. Similar trend has 
been demonstrated to occur, due to the absence of vis- 


cous effects, on a flat delta wing (Ref. 33). 

Finally, the inviscid results presented in figure 7 
also reveal the correlations between the surface pressure 
distributions obtained from structured (dash lines) and 
unstructured-grid (dash-dot lines) methodologies. The 
correlations clearly indicate an excellent agreement be- 
tween the results both in terms of magnitudes as well 
as the overall characteristics at all stations examined. 

Similar procedures were employed to compute the 
flow about the isolated MTVI fuselage at a higher an- 
gle of attack. The computational results based on both 
structured- and unstructured-grids are obtained at a = 
29.8°, Mqo ~ 0.4, Rft = 2.3 x 10 6 with available exper- 
imental data. The computed surface pressure distribu- 
tions obtained from these solutions are presented in fig- 
ure 8 along with the corresponding experimental data. 
Very similar conclusions can generally be drawn from 
these results as those identified previously for the flow 
conditions at a — 19.8°, M^ — 0.4, Rft — 2.4x 10 6 . In 
particular, note the excellent correlations between the 
structured- and unstructured-grid inviscid solutions. 

Forces and moment correlations: The predicted 
longitudinal aerodynamic characteristics are correlated 
with those measured experimentally over the isolated 
MTVI fuselage model in figure 9. The wind-tunnel data 
is presented over the entire range of flow-conditions 
tested, while the computational results are obtained 
for those conditions listed in table I. Similar to pre- 
vious figures (7, and 8), the results shown in figure 9 
are intended to provide a measure for assessing: 1)- 
the accuracy of the structured-grid viscous solutions in 
predicting the experimental data; 2)- the viscous effects 
from the structured-grid solutions; and 3)- the validity 
of the unstructured-grid inviscid solutions. 

The computed lift and drag coefficients obtained 
from the structured-grid viscous solutions (filled cir- 
cles) compares very well with the measured data (open 
circles), however the pitching moment appears to have 
been slightly overpredicted. The viscous effects on the 
total forces and moment can be identified by compar- 
ing the structured-grid viscous (filled circles) and in- 
viscid (filled squares) solutions. Despite the signifi- 
cant influence on the computed surface pressure dis- 
tributions (see figs. 7 and 8), the viscous effects ap- 
pear to cause only small change in the total lift, drag, 
and pitching moment coefficient. Finally, the longitudi- 
nal aerodynamic characteristics computed by the invis- 
cid unstructured-grid method (open squares) correlates 
reasonably well and those obtained from the structured- 
grid (filled squares) method, for the drag and pitching 
moment coefficient. However, the lift coefficients ap- 
pear to have been overpredicted by the unstructured- 
grid method. 
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Convergence & Performance characteristics 

The present computations are all performed on the 
Cray-YMP, located either at NASA Ames (Numerical 
Aerodynamic Simulation (NAS)) or at NASA Langley 
Research Center. The convergence history for both the 
viscous/inviscid structured- and inviscid unstructured- 
grid solutions, obtained at a = 19.8°, — 0.4, 

R/t = 2.4 x 10 6 , are shown in figure 10. The structured- 
grid viscous solution took about 2,600 cycles to reduce 
the residuals by 6 orders of magnitude and limit the 
Cl oscillations to a negligible level. The structured- 
grid inviscid solution took 1,500 cycles to reduce the 
residuals by about 8 orders of magnitude and limit the 
Cl oscillations to a comparable level as those of the 
viscous solution. The unstructured-grid inviscid solu- 
tion required only 500 cycles to reduce the residuals by 
about 4.5 orders of magnitude with a similar reduction 
in Cx, -oscillation characteristics. 

The effects of structured grid refinement on both 
the viscous and inviscid solutions were not investigated 
due to: 1)- the adopted viscous grid-strategy (designed 
to adequately resolve the geometrical complexity along 
with the corresponding flow physics) are comparable to 
those that have been found suitable (Ref. 27) for resolv- 
ing such vortical flows emanating from a sharp-edged 
configuration. 2)- the rational procedure employed for 
eliminating the grids in the boundary-layer region (de- 
termined by analyzing the viscous solutions from pre- 
vious step), to generate the inviscid grids. However, 
since the sensitivity of USM3D solutions to grid refine- 
ment has not been addressed for vortical flows, an ef- 
fort is made here to evaluate that effect for the present 
unstructured-grid analysis. 

Two new unstructured-grids are generated about 
the isolated MTVI fuselage; one coarser and the other 
finer with respect to the baseline grid. The coarse grid 
employed cell sizes that are globally about 15% larger 
while the fine grid used cell sizes that are about 10% 
smaller, compared to the baseline grid. Similar proce- 
dures are employed to obtained the solutions for the 
coarse and fine grids as those used with the baseline 
grid. The spanwise surface pressure distributions com- 
puted at a = 19.8°, = 0.4, are shown in figure 

1 1 for the same three fuselage stations used in previous 
correlations. The corresponding number of surface tri- 
angles as well as the total number of cells for each grid 
strategy are also denoted in the figure. The correlations 
indicate that the baseline grid provides adequate grid 
size and density to resolve the subject flow due to the 
small sensitivity of the corresponding surface pressure 
coefficients to further grid refinement. 

The computational efficiency and performance 
characteristics of the present unstructured-grid method 
are contrasted, in figure 12, against those of the 


structured-grid method obtained from both inviscid 
and viscous analysis. Particular attention should 
be paid to account for the disparity in the total 
number of cells used by each computational method 
(unstructured-grid has more than twice the number of 
cells) before interpreting the other statistics presented 
in the table such as the memory requirement and to- 
tal CPU time. As discussed earlier, this disparity in 
the total number of cells is primarily attributed to the 
constraint (i.e., mild cell stretching in the axial direc- 
tion) imposed by the present unstructured-grid genera- 
tion method. Improvements to this constraint, to allow 
for tetrahedral-mesh generation with a comparable ax- 
ial stretching as those of the structured-grids, should 
at least alleviate the existing differences (in terms of 
memory requirements and CPU-time) by a factor of 
about two. Despite the latter improvement however, 
the unstructured-grid technology is generally known 
to be more computationally intense, both in terms of 
memory requirements and the computational time, due 
to the random nature of the grid distribution. Nonethe- 
less, these computational resource requirements can be 
trade-offs for the quick turn-around time to generate 
meshes around complex configurations. 

Concluding Remarks 

Inviscid flow analyses are presented for vorti- 
cal flow prediction capability of an unstructured grid 
methodology. The accuracy of the results are ver- 
ified through a comparison with the computational 
predictions obtained with the structured-grid technol- 
ogy based on both the thin-layer Navier-Stokes with 
turbulent assumption as well as the inviscid formula- 
tions. The validity of the structured-grid turbulent so- 
lutions are first supported by a complementary set of 
wind-tunnel data. The viscous effects on vortical flow 
structures are then identified within the structured-grid 
framework. Subsequently, the structured-grid inviscid 
results are used to provide the bench-mark measures 
to evaluate the accuracy of the results obtained from 
the unstructured-grid methodology. The latter corre- 
lations revealed an excellent agreements between the 
computational results obtained from structured- and 
unstructured-grid methods in terms of off-surface flow 
characteristics and surface pressure distributions. The 
unstructured-grid inviscid solutions with the baseline 
grid were found to be insensitive to further grid refine- 
ment. The convergence history and performance char- 
acteristics obtained from each computational method 
are also contrasted which identifies their respective 
strengths and weaknesses in terms of resource require- 
ments. 
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(a) Complete configuration. 





lili 










(b) Isolated fuselage. 

Fig. L- MTVI wind- tunnel model, 
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Fig, 2.- Structured surface-grid for the isolated MTVI fuselage • 



a)- Farfield view 



Fig. 3.- Structured viscous (dotted lines) and inviscid (solid lines) flowfield grid strategies. 
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Fig. 4 - Unstructured surface-grid for the isolated MTVI fuselage • 




Fig. 5.- Unstructured inviscid flowfield-grid. 
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Contour Levels : min. =0,86, max. = 1.0, inc. =0.0175 


Fig. 6.- Crossflow normal! zed- total-pressure predictions based on structured and unstructured-grid methodologies; 

a= 19.8°, Moo = 0.4, R fi = 2,4x 10 6 . 
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Fig. 7.- Predicted surface pressure coefficients and correlations with data; a = 19.8°, — 0.4, Rjt — 2.4 x 10 6 



Local span, inches Local span, inches 


Fig. 8.- Predicted surface pressure coefficients and correlations with data; a = 29.8°, Moo = 0.4, Rf t — 2.3 x 10 6 
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Fig. 9.- Predicted longitudinal aerodynamic characteristics and correlations with data; M 0 0 = 0.4, Rfi ta 2 x 10 6 . 




Fig. 10.- Structured/unstructured-grid solutions convergence history; a — 19. 8®, M, oo = 0.4, Rjt — 2.4 x 10 6 . 
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# Surface A 

# Cells 

Coarse 

8,639 

361,122 

Baseline 

9,244 

567,424 

Fine 

13,350 

770,825 



Fig. 11- Unstructured-grid solution sensitivity to grid refinement; a — 19.8°, M = 0.4. 
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boundary faces 
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Memory requirements 
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Total CPU time 
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Fig. 12.- Structured/unstructured-grid methods performance characteristics; a = 19.8°, M ^ = 0.4, Rj t = 2.4 x 10 6 
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